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The incorporation of bismuth (Bi) in GaAs results in a large reduction of the band gap energy 
(E 9 ) accompanied with a large increase in the spin-orbit splitting energy (Aso), leading to the 
condition that Aso > Eg which is anticipated to reduce so-called CHSH Auger recombination 
losses whereby the energy and momentum of a recombining electron-hole pair is given to a second 
hole which is excited into the spin-orbit band. We theoretically investigate the electronic structure of 
experimentally grown GaBi^ Asi-^ samples on (100) GaAs substrates by directly comparing our data 
with room temperature photo-modulated reflectance (PR) measurements. Our atomistic theoretical 
calculations, in agreement with the PR measurements, confirm that E 9 is equal to Aso for x ~ 9%. 
We then theoretically probe the inhomogeneous broadening of the interband transition energies as a 
function of the alloy disorder. The broadening associated with spin-split-off transitions arises from 
conventional alloy effects, while the behaviour of the heavy-hole transitions can be well described 
using a valence band-anticrossing model. We show that for the samples containing 8.5% and 10.4% 
Bi the difficulty in identifying a clear light-hole-related transition energy from the measured PR data 
is due to the significant broadening of the host matrix light-hole states as a result of the presence of 
a large number of Bi resonant states in the same energy range and disorder in the alloy. We further 
provide quantitative estimates of the impact of supercell size and the assumed random distribution 
of Bi atoms on the interband transition energies in GaBi^Asi-^. Our calculations support a type-I 
band alignment at the GaBi^ Asi-^/GaAs interface, consistent with recent experimental findings. 



I. INTRODUCTION 

There is increasing interest in the highly mis- 
matched semiconductor alloy, GaBi-rAsi-^, both from 
a fundamental perspective^ and also because of its 
potential device applications!^! When a small frac- 
tion of As is replaced by Bi in GaAs, the energy gap 
Eg decreases rapidly, by « 90 meV when 1% of As 
is replaced by Bi. In addition, photo-reflectance mea- 
surements show that the energy separation, A$o, be- 
tween the spin-split-off valence band and the valence 
band edge also increases rapidly with Bi composition.^ 
This then leads to the situation in GaBi^Asi-^, where 
the spin-orbit-splitting energy can exceed the energy 
gap Aso > Eg for x > 9%. This is of significant po- 
tential benefit for telecomm lasers, because it could 
enable the suppression of the Auger recombination 
losses which dominate the thre shold characteristics of 
GalnAsP and AlGalnAs lasers PaiSI 

The strong band gap bowing in GaBi^ Asi-^ is sim- 
ilar to that observed in GaN^Asi-x, where E g de- 
creases by as much as 150 meV when 1% of As is re- 
placed by N. The band gap reduction in GaN^Asi-a; 
has been well explained using a band-anticrossing 
(BAC) modelP^ It is well established that replac- 
ing a single As atom by N introduces a resonant de- 
fect l evel a bove the conduction band edge (CBE) in 
GaAs EzUHl This occurs because N is considerably more 
electronegative and is also 40% smaller than As. The 
interaction between the resonant N states and the 
CBE of the host GaAs matrix then accounts for the 
observed rapid reduction in E g . It was proposed^, 
and later confirmed theoretically^^ that replacing As 



by Bi can introduce a similar BAC effect. Like N 
and As, Bi is also a group- V element, which lies be- 
low Sb in the periodic table. Bi is much larger than 
As, and is also less electronegative. It should there- 
fore be expected that any Bi-related resonant defect 
levels should lie in the valence band (VB) and that, 
if an anticrossing interaction occurs, it will occur be- 
tween the Bi-related defect levels and the valence band 
edge (VBE) of the GaAs matrix. Recent sp 3 s* tight- 
binding calculations which we have performed confirm 
the existence of such a band-anticrossing interaction, 
with the calculated variation of E g and of A go in ex- 
cellent agreement with experiment P 

Although the BAC model provides an accurate de- 
scription of some of the key features in the band struc- 
ture of highly mismatched alloys, it includes several 
critical simplifying assumptions. It assumes for in- 
stance in GaBLuAs!^ that each Bi atom introduces 
resonant defect levels with constant energy Esi, and 
then treats the effect on the band structure in terms 
of the interaction between these states and the host 
matrix valence band. The BAC model therefore effec- 
tively assumes that all Bi atoms are in an identical en- 
vironment. Although this may be approximately true 
for very low Bi compositions, where most Bi atoms are 
widely separated from each other, this assumption can 
be expected to break down with increasing Bi compo- 
sition. Firstly, there should be an increasing number 
of Bi pairs and clusters formed with increasing Bi com- 
position, where a Bi pair consists of a Ga atom that 
has two Bi neighbours, with clusters then containing 
larger numbers of contiguous Bi atoms. These pairs 
and clusters introduce defect levels which lie above the 
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isolated Bi state energy, Esi and which also interact 
with the valence band states 21 . In addition, isolated Bi 
atoms also see an increasingly disordered local envi- 
ronment, leading to an inhomogeneous broadening of 
the associated resonant defect level energies. We in- 
vestigate here the effects of these two different types of 
disorder on the valence band structure of GaBi-rAsi-^ 
grown pseudomorphically on GaAs, using the sp 3 s* 
tight-binding model introduced in Ref. 2. 

This paper presents a combined theoretical and ex- 
perimental analysis of the effects of disorder on the 
electronic structure of GaBi^Asi-a; epilayers grown 
under compressive strain on (100) GaAs substrates. 
Such layers are of interest for a number of reasons. 
Firstly, many applications such as semiconductor 
lasers and related devices require the growth of such 
layers. Secondly, the growth of these strained struc- 
tures has an interesting impact on the band structure, 
splitting the degeneracy of the heavy-hole and light- 
hole states at the valence band maximum, with the 
splitting between the two sets of states initially in- 
creasing for instance at a rate of about « 75 meV per 
% lattice mismatch for a compressively strained In- 
GaAs layer grown pseudomorphically on GaAs. We 
note that the compressive strain in a GaBi x Asi_ x 
pseudomorphic layer should also split the four-fold de- 
generacy of the Bi resonant states, with the two states 
with heavy-hole symmetry shifting upwards in energy, 
while the light-hole states shift downwards. With in- 
creasing strain, the GaAs host matrix light-hole states 
should start to pass through the Bi heavy-hole- like res- 
onant states. The intrinsic disorder in GaBiAs then 
lead to mixing between these two types of state. 

The GaBi^Asi-j; samples studied were grown by 
molecular beam epitaxy with Bi concent rations of 
x = 2.3%, 4.5%, 8.5%, and 10.4%PMl2l The sam- 
ples are fully strained (pseudomorphic) to the GaAs 
substrate, as confirmed by x-ray diffraction data in- 
dicating the GaAs in-plane lattice constant for the 
GaBi^Asi-z layersP^HIS The theoretical analysis was 
undertaken using the sp 3 s* tight-binding Hamilto- 
nian which we have developed for GaBiAs^ to investi- 
gate the electronic structure of large randomly disor- 
dered supercells. The electronic structure is studied 
experimentally by photo-modulated reflectance (PR) 
spectroscopy^, which is considered to be an excellent 
technique due to its sensitivity to critical point tran- 
sitions in the band structure.^ Further details of the 
experimental procedure are described in Batool e t al® 
while the sample details can be found in Lu et al! 16 * 17 ! 

Three sets of features are found in the PR spectra, 
associated with transitions between the conduction 
band minimum at T and valence states which respec- 
tively include host matrix heavy- hole (HH) , light-hole 
(LH) and spin-split-off (SO) character. Our calcula- 
tions provide a detailed understanding of the broad- 
ening observed in the PR measurements for each of 
these features due to the presence not just of isolated 
Bi atoms but also of pair and cluster states. We find 
different factors affecting each of the three transitions, 
with the calculated data being in generally good agree- 
ment with the measured spectra. In order to test the 
effects of disorder, we undertake a set of supercell cal- 



culations for each of the 4 samples investigated, where 
the supercell size is varied from 1000 to 8000 atoms, 
with several calculations with different random distri- 
butions of Bi atoms being undertaken for each super- 
cell size. We find for each supercell size that the cal- 
culated value of the energy gap separating the lowest 
conduction and highest (HH-like) valence states varies 
for a given Bi composition x depending on the actual 
distribution of Bi pair and cluster states in the su- 
percell considered, consistent with the inhomogeneous 
broadening of the lowest energy transition in the PR 
spectra. The calculated transition energy to the SO 
band is largely independent of the supercell size and 
configuration, but we find that the GaAs host matrix 
SO state character is spread over a number of supercell 
states, contributing also to an inhomogeneous broad- 
ening of the transition energy, which increases with 
increasing Bi composition x. Finally we consider tran- 
sitions to states which include GaAs LH zone-centre 
character. Because these states become degenerate 
with a broad distribution of Bi-related resonant states 
in a strained layer, we find that transitions to the LH 
states become very strongly perturbed with increasing 
Bi composition, with the GaAs LH T character distri- 
bution varying both with supercell size and with the 
particular random distribution of Bi compositions. At 
higher Bi compositions, no single feature can be asso- 
ciated with the GaAs LH states, consistent with the 
difficulty in fully fitting the LH part of the spectrum 
at higher Bi compositions. 

The remainder of the paper is organised as follows. 
We give a brief overview of the theoretical and ex- 
perimental methodology in section |H} The main re- 
sults are then presented in section |III| We begin by 
presenting a general comparison between the exper- 
imental and theoretical results in section MI Al This 
is followed in section III B| by a more detailed anal- 
ysis of the effect of supercell size and of assumed Bi 
atom distribution on the calculated spectra. We con- 
sider supercell sizes of 1000, 4096, and 8000 atoms for 
each of the four grown samples. We find for a fixed 
Bi composition that a larger supercell size results in 
an enhanced inhomogeneous broadening in particular 
of the LH character. The results clearly highlight the 
importance of performing large supercell calculations 
to model realistic dilute alloy environments. In sec- 
tion |III C[ we quantify the impact of the random dis- 
tribution of Bi atoms on the broadening of transition 
energies by performing 8000 atom supercell calcula- 
tions with four different random distributions of Bi 
atoms for each of the four compositions considered. 
Our results indicate a significant broadening of the 
transition energies due to the random distribution of 
atoms, showing that the local band edge energies in a 
GaBiAs alloy are very sensitive to the Bi distribution 
in the given region. In section |III D we summarise 
the key physical insights gained by comparison of our 
atomistic simulations with the measured PR spectra. 
In section [ill E| we consider the band edge alignment 
in compressively strained GaBi^Asi-^/GaAs samples. 
The band alignment at the GaBi^Asi-^/GaAs inter- 
face remains uncertain in the literature with reports 
suggesting that the conduction band offset is of type 
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Ij 2 ' 20 ' 21 ' or type IlfSl or nearly flatP Our calculations 
support a type-I band alignment for all the samples 
considered. The overall conclusions of our work are 
summarised in section ITVl 



II. METHODOLOGY 

The atomistic calculations are performed using a 
tight-binding (TB) model, where the alloy atom posi- 
tions are determined using a valence force field (VFF) 
strain energy minimization scheme and the electronic 
structure is then calculated using an sp 3 s* TB Hamil- 
tonian, including bond-length and bond-angle depen- 
dent interaction parameters. The VFF and TB mod- 
els are implemented inNanoElectronic MOdeling Sim- 
ulator (NEMO 3-D)P124| T h e values of the tight- 
binding parameters for GaAs and GaBi, as well as the 
values of the force constants (a and /3) for the VFF 
model are reported in Ref. 2. The VFF parameters for 
GaBi were obtained by fitting to first principle calcula- 
tions of GaBi elastic constants.^ We further verifiecP 
that the relaxed Ga-Bi bond lengths in GaBiAs cal- 
culated using the VFF model were in good agreement 
with the results of the x-ray absorption spectroscopy 
measurements thereby confirming the validity of 
the VFF mdoel employed. The TB and VFF models 
have previously been applied to study free-standing 
GaBi^Asi-z and GaBi^Pi-a; supercells containing up 
to 8000 atoms. These calculations showed that Bi- 
related resonant states interact with GaAs host ma- 
trix valence band states in GaBi^Asi-^ via a composi- 
tion dependent band anti-crossing (BAC) interaction, 
which we calculate to shift the highest valence band 
edge upwards in energy by « 55 meV per % of Bi 
replacing As for low Bi compositions. This BAC in- 
teraction is analogous to the BAC interaction between 
nitrogen (N) defect states and GaAs conduction band 
states in GaN^Asi-a; alloys.^ The calculated compo- 
sitional dependence of the band gap and spin-orbit- 
splitting energies demonstrated excellent agreement 
with the available experimental data sets including 
Bi samples with compositions up to 10%. The ex- 
perimentally observed large bowing of the band gap 
energy was explained in terms of two factors: (i) a 
linear decrease of ss 30 meV per % Bi in the lowest 
conduction band edge energy (E c i) which occurs due 
to a conventional alloy shift and (ii) a large non-linear 
increase in the highest valence band edge energy (Ei,i) 
of w 30-55 meV per % Bi due to the band anti-crossing 
interaction. The energy of the spin-split-off band (SO) 
was calculated to shift downwards in energy slightly, 
exhibiting a linear decrease of w 5 meV per % Bi. It 
should be noted that no evidence was found for any 
BAC interaction involving the SO band, as earlier as- 
sumed by Alberi et al^ 

In this paper, we apply our TB model to investigate 
compressively strained GaBi^Asi-^ samples, grown 
epitaxially on a (100) GaAs substrate. In the VFF 
relaxation scheme, we keep the in-plane lattice con- 
stant for all supercells considered fixed to that of the 
GaAs lattice constant, while allowing the supercell to 
relax freely along the growth direction. We use peri- 



odic boundary conditions in all three spatial directions 
for the calculations of electronic structure. For each 
calculation undertaken, the Bi atoms are placed on As 
sites, selected randomly in the supercells. 

In order to compare the theoretical calculations 
with the experimental data, we calculate the distri- 
bution of the host matrix T-character, Gr(E), across 
the alloy valence band states. The calculation of the 
T-character is carried out by projection of the six high- 
est host matrix zone centre valence band states (two 
HH, two LH and two SO) onto the full spectrum of 
the energy states in the alloy (GaBi^Asi-z) supercell 
under consideration. In what follows, we use the sub- 
scripts /, and k, 1 to index the unperturbed host and 
Bi-containing alloy states, respectively. The Gr(E) 
spectrum is then calculated by projecting the unper- 
turbed GaMAsM supercell states at T, \ipi,o), onto 
the spectrum of the GaMBiLAsA/_L alloy supercell, 
{E k ,\1>k,i)Y- 



g(Ei) 

G r(£) = E E fr,uHEk-E) (1) 
fc i=i 

fr,H = l(^,il^,o)| 2 (2) 

Here 6 (Ek — E) is the Dirac delta function centred 
on E^, g{Ei) is the degeneracy of the host matrix 
band having energy Ei at the T point in the Brillouin 
zone and fr,kl is the fractional T character of the 
alloy state \ipk,i)- Since the PR measurements are 
performed at room temperature (300 K) and the TB 
calculations are parametrised to the low temperature 
band structure, we shift the energy scale of our theo- 
retically calculated T-character plots at K by 0.098 
eV to convert them to room temperature values. 

The analysis of Gr(E) plots has provided a very 
useful approach to understand t he ch aracter and evo- 
lution of alloy band edge statesP^H Analysis of such 
plots confirms the presence of the band anti-crossing 
interaction in ordered supercell calculations of dilute 
bismide alloys 

Turning to the experimental data, the measured 
PR spectral line shapes were least-squares fitted using 
Aspnes third derivative functionals of the forrrP^ 



=Re[Ce w (E-E g +t 1 y n } (3) 

where C and 9 are amplitude and phase variables, E 
the energy of the probe beam, E g the critical point 
transition energy, and 7 a broadening parameter. We 
investigated the effect of choosing several of the sug- 
gested values for the line shape exponent factor np^i 
but found that this had little influence on the fitted 
transition energies. Here, we present the results for 
fits using n = 3. In the next section |III A| we use 
Eq. [3] to fit the measured PR data and identify the 
interband transition energies. 
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TABLE I: The values of the broadening parameters 7 used in Eq. [3] to fit the PR spectra plotted in Figs. [T]and[2] 



Bi composition 


JHH 


"JLH 


ISO 


(%) 


(meV) 


(meV) 


(meV) 


2.3 


82.1 ± 5.3 


51.9 ± 3.5 


30.9 ± 0.3 


4.5 


31.3 ± 1.2 


31.0 ± 2.0 


22.2 ± 1.3 


8.5 


31.4 ± 0.5 


82.1 ± 7.0 


31.9 ± 2.2 


10.4 


31.4 ± 1.0 


31.4 ± 1.0* 


38.4 ± 3.2 



*Here j LH was constrained to equal ^hh in the 
fit for 10.4% Bi, due to an extra feature in the 
LH energy range which is difficult to model. 



III. RESULTS AND DISCUSSIONS 
A. Comparison of Theory and Experiment 

Fig. [I] shows the measured PR spectra (black dots) 
near to and just above the energy gap for the four 
GaBi^Asi-z samples considered, with compositions x 
= 2.3%, 4.5%, 8.5%, and 10.4%. The measured PR 
data have been fitted using Eq. [3j with the fits shown 
in the figure by the solid red lines. The values of 
the broadening parameter 7 used in Eq. [3] to obtain 
best fits to the experimental data are listed in table 
|TJ These values are directly related both to inhomo- 
geneous broadening of the transition energies due to 
alloy disorder effects in the grown samples, as well 
as to any composition fluctuation effects. The fitted 
7 values will be compared to the broadening of our 
theoretically calculated T-character plots later in this 
paper. 

For each composition case shown in Fig. [T] cumu- 
lative HH+LH r character spectra are also plotted, 
using the same energy scale (green bars) and where 
each peak in the theoretical spectra is associated with 
a transition between the lowest conduction state and 
a given valence state. The T-character plots are com- 
puted from 8000 atom supercell calculations, with the 
Bi atoms randomly distributed within each supercell 
considered. Here we consider only one random distri- 
bution of Bi atoms to compare against experiment for 
each composition. The results here therefore highlight 
alloy disorder effects such as mixing of alloy valence 
band states with Bi pair and cluster related states, 
as a function of the Bi composition. Later in section 
|f II C[ we will investigate the second important factor 
that contributes in the inhomogeneous broadening of 
the transition energies by performing calculations for 
four different random distributions of Bi atoms for 
each grown sample. These later results will show that 
the specific distribution of Bi atoms within a super- 
cell can further alter the calculated energy spectra by 
introducing different Bi local environments. This can 
result in a redistribution of T-character and shifts in 
the transition energies, thus providing a further con- 
tribution to the inhomogeneous broadening of the cal- 
culated spectra. 

We find that the alloy conduction band (CB) V- 
character is largely concentrated on the lowest energy 
conduction state of the alloy, which retains more than 
95% r character for all of the supercells and Bi compo- 
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FIG. 1: Room temperature photo-reflectance (PR) spectra 
for the four compressively-strained GaBi^Asi-^ samples, 
in the region of the fundamental band gap highlighting 



the E^" 1 = E i- 



Eui and E^"" 2 = E c i— E„2 transitions. 



The black circles are from measured values and the red 
lines are the fits obtained by equation [3] The green bars 
are the calculated GaAs HH+LH T-character plots for the 
alloy valence band states, computed from 8000 atom su- 
percell calculations. The values of the HH and LH related 
transition energies are selected from the F character plots 
as the energies with the highest HH and LH F character 
projections, respectively. The dashed lines are a guide to 
the eye and indicate the PR fitting transition energies. 



sitions considered. This is consistent with our previ- 
ous free-standing GaBi^Asi-a; supercell calculations, 
where the conduction band T-character was found to 
be more than 99.9% concentrated on the lowest energy 
stateP This confirms that we can use a conventional 
alloy model to describe the shift in the conduction 
band edge energy with alloy composition. Therefore, 
we do not show plots of the CB T-character in this 
paper, and use the calculated lowest conduction band 
edge energy, E c i , to compute the interband transition 
energies. 

The two composition dependent features in the PR 
plots of Fig. [I] are assigned for the 2.3% and 4.5% 
samples to the E^" 1 = E ci - E„i and E^" 2 = E ci - 
E„2 transitions. This assignment is supported by the 
T-character plots, where the lower energy feature is 
associated with transitions to HH-like valence states 
(E v i) and the higher energy feature with transitions 
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FIG. 2: Room temperature photo-reflectance (PR) spectra 
for the four compressively strained GaBi^ Asi-a; samples 
in the region of the spin-orbit split off transition (E^ so 
= Ed— Eso)- The black dots are from measured val- 
ues and the red lines are the fits described by equation [3] 
The green bars are the calculated GaAs SO F-character 
plots for the alloy valence band states, computed from 
8000 atom supercell calculations. The upper arrows in 
each panel indicate the PR fitting energies and the lower 
arrows mark the calculated energies from the F character 
plots. The prominent feature around 1.77 eV in the PR 
spectra is due to the SO related transition from the GaAs 
substrate. 



to a state {E V 2) with significant LH character. The 
fittings obtained from Eq. [3] for these samples, shown 
by the red lines are in good agreement both with the 
theory and with experiment. 

As the Bi composition is increased to 8.5% and 
10.4%, increasing disorder, including the presence of 
a large number of Bi pair and cluster states, results in 
a large inhomogeneous broadening of the T-character 
plots over the higher valence band states. This makes 
the E^ 1 ' 2 transition difficult to identify in the PR spec- 
tra, as can be observed from the poor agreement be- 
tween the experimental data (black circles) and the 
fitted values of ^ (red lines) for these two com- 
positions. The corresponding broadening of the T- 
character plots is highlighted by the dashed (blue) 
ovals in Fig. [l] consistent with the difficulty in fitting 
the high energy feature in the measured PR spectra. 
It should be noted that the low energy transition in- 
volving E c x and E„i does not require an equivalent 
broadening either in theory or experiment. This is 
due to the fact that the energy of the E„i state shifts 
upwards in compressively grown samples, and for high 
Bi compositions (8.5% and 10.4%), it can then still re- 
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FIG. 3: Plot of the room temperature transition ener- 
gies, Ef" 1 , Ef" 2 , and Ef so , obtained from fitting the PR 
spectra (black solid lines) and calculated from the atom- 
istic simulations (broken red lines). The values are directly 
extracted from Figs. [I] and [2] 



tain significant host matrix HH T character. 

Fig. [2] shows the measured PR spectra (black cir- 
cles) , fits of Eq. [3] (red lines) , and the calculated GaAs 
SO band T-character plots (green bars) for the spin- 
split-off transition energy (E^ so = E cl - E so ). The T- 
character for the SO band is spread over several states, 
which are closely spaced in energy. This makes it dif- 
ficult to choose a particular energy value for the SO 
transition energy. We select the SO transition energy, 
illustrated by the lower vertical arrow in each panel, 
by calculating the weighted average of the energies 
with their corresponding T-character. The transition 
energies obtained from fits to the experimental data 
using Eq. [3] are shown by the upper vertical arrows in 
each case. 

Fig. [3] compares the calculated values of the 
three composition dependent transition energies E^" 1 , 
E^" 2 , and E^ so , with the fitted PR transition ener- 
gies extracted from Figs. [1] and |J For the four Bi 
compositions considered here, our theoretical values 
of the transition energies are overall in good agree- 
ment with the experimental values. The fittings of 
the E^" 2 and E^ so energies for the 10.4% sample are 
less good, due to a significant broadening of the LH 
and SO T-character which makes it difficult to pick a 
particular energy accurately from the PR data. 

From the results of Fig. [3J we can further compute 
the spin-orbit-splitting energy Aso = E„i — E50 as 
A so = E^ so - E^" 1 . Fig. S| shows the plots of E g = 
E^" 1 and of Aso obtained both from our TB calcu- 
lations and from the experimental PR measurements. 
The two results are in good agreement and the cal- 
culated value of the Bi fraction where E ff = Aso (~ 
9.6%) is in excellent agreement with the value (« 9.0 
± 0.2%) obtained from the fits to the PR data. This 
crossing is technologically important for the design 
and realisation of photonic devices, due to the possi- 
bility of suppressing at higher bismuth concentrations 
the CHSH Auger recombination loss process, which 
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FIG. 4: Comparison of the experimental and theoretical 
values of the energy gap E 9 and spin-orbit-splitting energy 
Aso obtained from the results of Fig. [3] The values of the 
Bi composition where E 9 = Aso obtained from theory and 
experiment are also shown. 



involves hole excitation from the highest valence band 
into the spin-split-off bandP^LOI] 



B. Effect of Supercell Size 

The size of supercell chosen in any theoretical study 
of dilute impurity alloys is an important calculational 
consideration, because it can significantly impact the 
accuracy of the results. By performing ordered super- 
cell calculations, it has previously been shown that a 
minimum supercell size of 4096 atom s is required to 
realize the dilute impurity limfPESl in GaBi^Asi-^, 
GaBi^Pi-a;, and InBi^Asi-^ alloys. The larger su- 
percell size also increases the accuracy of results by 
suppressing spurious cell-to-cell interactions between 
impurities that may arise due to the periodic bound- 
ary conditions.^ 

In this section, we investigate as a function of the 
supercell size the impact of alloy disorder on the 
broadening of T-character plots and on the values of 
the transition energies. We consider three supercell 
sizes containing 1000, 4096, and 8000 atoms respec- 
tively and compare the HH+LH T-character plots for 
all four Bi compositions in Fig. [5] 

For the 1000 atom supercells, we calculate two dis- 
tinct T-character peaks associated with the two high- 
est valence band states for all four Bi compositions 
considered. This is clearly not consistent with the 
measured PR data in Fig. [l] where for the 8.5% and 
10.4% samples, the higher energy transition corre- 
sponding to the E„2 state is significantly broadened 
and becomes hard to define at a single energy value. 

The larger 4096 and 8000 atom supercell calcula- 
tions are both able to capture this effect where only 
the highest valence band state E^i is distinct in the 
T-character plots and the rest of the T-character, in 
particular the TH-related character, is spread over a 
large number of valence states, making it difficult to 



associate a particular energy value with E„2- This is 
because for the 1000 atom supercells, the lower num- 
ber of Bi atoms results in the formation of fewer pair 
and cluster states, and the density of states is also 
lower due to the presence of fewer folded bands near 
to the GaAs valence band maximum at the T pointP 
Therefore, the interband transition involving E.„2 in 
theory remains relatively sharp. From these results, 
we conclude that large supercell calculations are cru- 
cial to accurately mimic a dilute impurity alloy. 

Another discrepancy in the 1000 atom calculations 
is prominent from the T character plots for the 8.5% 
and the 10.4% compositions in Fig. [5j where the 1000 
atom supercell calculations result for both composi- 
tions in the E„i state being at significantly higher en- 
ergy than in the 4096 and 8000 atom calculations. 
This indicates that the small size of the supercell can 
overestimate the reduction in the band gap energy 
at large Bi compositions. Further calculations which 
we have undertaken demonstrate that there can be 
a strong variation in the E v i energy obtained from 
different 1000 atom supercell calculations, where sta- 
tistical variations in the number of pairs or clusters in 
the supercell can have a significant impact on the cal- 
culated energy spectrum. We will see below that the 
calculated energy of the highest valence state does not 
vary significantly between different 8000 atom super- 
cell calculations. Overall, changing the supercell size 
from 4096 atoms to 8000 atoms provides a small re- 
duction in the variation between different T-character 
plots and therefore a choice can be made between 
these two supercell sizes, based on the trade-off be- 
tween desired accuracy and computational time. 

C. Effect of Random Distribution of Bi Atoms 

Another important factor that can significantly in- 
fluence the theoretically calculated electronic struc- 
ture of dilute impurity alloys is the particular random 
distribution of impurity atoms included in the simu- 
lated supercells. Different random distributions of Bi 
atoms will give varying numbers and types of Bi pairs 
and clusters in the supercell. This can in turn for 
smaller supercells introduce considerable variations in 
the T-character plots and in the values of the inter- 
band transition energies. 

In this section, we perform a series of 8000 atom 
supercell calculations for each of the four Bi composi- 
tions to investigate the effect of statistical variations 
on the electronic properties of the alloy. For each com- 
position, we consider four different random distribu- 
tions (RDs) of Bi atoms labelled as RD1, RD2, RD3, 
and RD4, which are obtained by varying the seed val- 
ues of a random number generator that determines 
the nature of anion (either Bi or As) at a given atomic 
location in the supercell. By using different random 
number seed values, we ensure that all of the four 
supercells (RD1-RD4) have different arrangements of 
the Bi atoms, resulting in different numbers and types 
of Bi pairs and clusters. 

Fig.|6]plots the sum of the HH and LH T-character, 
Gr {E) , for the four random distributions of Bi atoms 
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FIG. 5: For each Bi composition of the experimentally grown samples, we calculate and plot combined HH and LH F 
character for three alloy supercells consisting of 1000, 4096, and 8000 atoms. For the same Bi composition, a larger alloy 
supercell will have more Bi atoms and therefore an increased number of pair and cluster configurations. This is reflected 
by the wider F character distribution for larger supercells. The significant overestimation of the alloy VBE energy in 
1000 atom supercell calculations for 8.5% and 10.4% samples is due to artificially enhanced cell-to-cell interaction of pair 
and cluster states in the structures considered. 
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FIG. 6: For each of the four Bi compositions, we consider four random distributions of Bi atoms in the supercell, labelled 
as RD1, RD2, RD3, and RD4. Each random distribution introduces different pair and cluster states and therefore 
exhibits different broadening of the combined HH and LH T-character. The calculated data are shown for 8000 atom 
supercell calculations. We also specify the values of individual HH and LH contributions in the few highest energy bars 
of the plot as (HH,LH). If (,LH) or (HH,) is mentioned for a particular bar, that indicates that HH or LH contribution 
in that bar is less than 0.1%. The zero of energy is taken in all cases at the GaAs valence band maximum. 
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considered at each composition. We specify the sepa- 
rate contributions from the GaAs HH and LH states 
in the label (HH,LH) for the highest few states, in or- 
der to highlight the distribution of the HH and LH 
character over the alloy states. It should be noted 
that (,LH) and (HH,) refers to cases where a particu- 
lar alloy state has less than 0.1% HH or LH character, 
respectively. 

The impact of variations between different random 
distributions is more prominent for the two lower Bi 
compositions, 2.3% and 4.5%. This can be mainly at- 
tributed to two reasons: (i) For the lower Bi composi- 
tions, there is a significantly smaller number of pairs 
and clusters in a supercell, and therefore the impact 
of any variation in their type or number is more pro- 
nounced. Similar conclusions were presented earlier^ 
based on free-standing supercell calculations, where 
the addition of a single pair or a single Ga-centred 
3 Bi cluster in a 4096 atom Ga204sBiA/As204s-A/ su- 
percell (M=2, 3) resulted in very large variations in 
the valence band edge energies, (ii) The highest va- 
lence band edge energy is very close to the energies of 
the pair/cluster states and the stronger state mixing 
which this allows therefore results in the broadened T- 
character distributions. This behaviour is confirmed 
by close examination of RD1 for 2.3% Bi and RD1 and 
RD3 for 4.6% Bi. These are the distributions which 
give rise to the highest valence band maximum ener- 
gies, and also have the lowest F character associated 
with the highest state, reflecting the contribution of 
(higher energy) pair or cluster states in each of these 
cases. The large ^hh value required to fit the PR 
spectrum of the 2.3% sample may then include a con- 
tribution due to mixing of the host T states with pair 
and cluster states which lie close to the band edge at 
this composition. 

For the larger Bi compositions, 8.5% and 10.4%, the 
highest valence band edge energy increases well above 
the energies of the pair/cluster states and therefore 
a sharp peak is observed for E v \ in the T-character 
plots. This is quite evident for the 10.4% composition 
case, where 55%, 59%, 52%, and 61% HH character is 
found on the highest valence band state for the four 
random distributions. For these compositions, the LH 
T-character is in each case distributed over a broad 
range of valence band states, so that it then proves 
very difficult in the 10.4% cases to identify any single 
state as the LH-like valence band maximum, consis- 
tent with the difficulty in fitting the PR spectra in 
this energy range. 



D. 



Summary of the Alloy Disorder Effects and 
Comparison with the PR fits 



In the previous two subsections III B and III C 



we 

have investigated the calculated impact of alloy dis- 
order when varying two important parameters: su- 
percell size and the assumed (random) distribution 
of impurity atoms. Both of these factors have been 
found to significantly impact the calculated electronic 
structure of GaBi^Asi-j, alloys, and therefore require 
careful attention. Below, we summarize our results 



and the key physical insights gained in the previous 
subsections concerning the electronic properties of the 
alloy: 

• A very small supercell size can result in arti- 
ficially enhanced alloy disorder effects on the 
band edge energies, and can therefore introduce 
significant errors in the calculated band gap en- 
ergy. Small supercells also fail to correctly cap- 
ture the broadening of LH character observed 
in experimental PR data for higher Bi compo- 
sitions. Our calculations show that a supercell 
size of 4096 or larger is required to accurately 
model dilute bismide alloys. 

• The effect of atomistic randomness remains im- 
portant at low Bi compositions (2.3% and 4.5%) 
even in larger supercells, where variations in the 
local Bi environments can significantly impact 
the calculated distribution of T-character, con- 
sistent with the large 'Jhh value required to fit 
the 2.3% PR spectrum. 

• At higher Bi compositions (8.5% and 10.4%), 
the effect of atomistic randomness is less marked 
at the valence band maximum, because the E„! 
energy is shifted well above the pair/cluster 
state energies. 

• Overall, the broadening of the theoretically cal- 
culated T-character plots correlates well with 
the values of the broadening parameter 7 (pro- 
vided in table [I]) used to fit the PR spectra 
via Eq. [3] The large values of the 7 parame- 
ter used to fit the PR data for the 2.3% sam- 
ple may include a contribution due to compo- 
sition variations in the grown sample, as well 
as a contribution due to state mixing with pair 
and cluster states. For the three higher compo- 
sition samples, a constant value of ^hh can be 
chosen to fit the HH-related transition. This 
is consistent with the theoretically computed 
T-character plots, where the HH-character is 
largely present on the alloy VBE state and does 
not show much spreading even for different ran- 
dom distributions of the Bi atoms, as evident 
from the plots in Fig. [6j 

• The value of the jso broadening parameter in 
table [T] for the SO band slightly increases with 
Bi composition. This is consistent with the T- 
character plots shown in Fig. [2j which show an 
increase in broadening of the T-character with 
increasing composition. In order to quantify this 
effect, we estimate the calculated broadening of 
the SO r character by calculating the energy 
range AE around E50 which includes 70% of 
the SO character. The values obtained for 2.3%, 
4.5%, 8.5%, and 10.4% Bi compositions are <1 
meV, 2.7 meV, 5.7 meV, and 13.8 meV, clearly 
indicating an increasing broadening of SO F- 
character with increasing Bi composition, con- 
sistent with conventional alloy disorder effects. 

• The behaviour of the ^lh parameter values for 
the LH-related transition energy is also well ex- 



!) 



plained by the theoretical analysis. Table |T] 
shows a significant increase in jlh for the 8.6% 
Bi sample, while a good fit could not be obtained 
to the experimental data for the 10.4% sample, 
due to the presence of extra features in the mea- 
sured data. This is consistent with our TB cal- 
culations that show a significant broadening of 
the LH T-character over a series of alloy states 
at higher compositions, making it very difficult 
to reliably choose a particular energy value for 
the E„2-related transition energy in the 10.4% 
sample. 



E. Conduction Band Edge Deformation 
Potential and Type-I Band Edge Alignments 

Fig. [7] compares the Bi composition dependent 
band edge energies for the lowest conduction band 
edge (E c i) and the highest two valence band edges 
(E„i and E„ 2 ) for the compressively strained and un- 
strained (free-standing) GaBi x Asi_ x supercells, cal- 
culated from our atomistic tight binding model. The 
compressive hydrostatic component of strain pushes 
the conduction band edge energy upwards and the bi- 
axial strain leads to an increasing splitting between 
the two highest valence band states with increasing x. 

We calculate a strain induced net decrease in E cl of 
13 meV per % Bi. Using 5E C — a c e# where 5E C is the 
change in the conduction band edge energy and e# 
is the hydrostatic strain, we can estimate the value 
of CB deformation potential a c for the GaBi 2 ,As 1 _ a , 
alloy. Using the values from Fig. [7] at 10.4% Bi com- 
position, 8E C = 0.12 eV. The hydrostatic strain is en 
w -0.013. Therefore the calculated value of a c for the 
compressively strained GaBi^Asi-^ is -9.4 eV, which 
is comparable to, but of larger magnitude than the val- 
ues -7.17 eV and -5.08 eV reported earlier for GaAs 
and InAs materials.^ 

The determination of the alignment of conduction 
and valence bands is critical for the design and opti- 
mization of optoelectronic devices. The band align- 
ment at the GaBi x As 1 _ a; /GaAs interface remai ns un- 
certain in the literature with reports of type Ij 2 l 2 ° l 21 l 
type IIpl and nearly fla1p2l conduction band offsets. 
We calculate here that the lowest conduction band 
edge energy moves up by ps 13 meV per % Bi with 
increasing Bi composition. This upward shift is not 
large enough to overcome the « 28 meV per % Bi 
decrease in energy due to the conventional alloy-type 
shift. Therefore for all the bismuth compositions con- 
sidered, the conduction band edge of GaBLjAsi-a; re- 
mains below the conduction band edge of the GaAs 
barrier material. The valence band edge, on the other 
hand, increases in energy due to the BAC interac- 
tion effect in the free-standing supercells and com- 
pressive strain further enhances this upward shift, so 
the top most valence band edge of the GaBi x As 1 _ a; 
alloy is always well above the valence band edge in 
the GaAs barrier material. At 10.4% Bi composition 
for which A so > E g , we calculate band offsets of 
160 meV and 400 meV for the conduction and va- 
lence band edges with respect to the GaAs bulk band 
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FIG. 7: The Bi composition dependent shifts in the band 
edges energies E c i, E„i, and E„2 of the GaBi^Asi-^ al- 
loy are shown, obtained from the compressively strained 
and free-standing 8000 atom supercell calculations at low 
temperature. The strain shifts the lowest conduction band 
edge to higher energy, reducing its slope from 28 meV/%Bi 
to 15 meV/%Bi. The splitting between the two highest 
valence band edges is also increased due to compressive 
biaxial strain. 



edges, respectively. From this discussion, we con- 
clude that GaBi^Asi-j; alloys grown pseudomorphi- 
cally on a GaAs substrate have a type I band align- 
ment at the GaAs barrier interface, i n acc ordance 
with more recent experimental findings |2llS! The in- 
crease in the calculated value of a c for GaBi^Asi-a; 
and the calculated type I band alignment are, in par- 
ticular, consistent with the results of hydrostatic pres- 
sure measurements undertaken on GaBi^Asi-^/GaAs 
LED structures^ These measurements found that the 
electro-luminescence intensity from GaBiAs decreased 
with increasing pressure, accompanied by an increase 
in the luminescence from the GaAs layer, consistent 
with increased carrier overflow with increasing pres- 
sure, owing to a type I offset and a larger conduction 
band deformation potential in the GaBiAs layer, than 
in the surrounding GaAs layers. 



IV. CONCLUSIONS 

In summary, we have investigated alloy disorder 
effects on the electronic structure of compressively 
strained GaBi^Asi^a; samples grown on a GaAs sub- 
strate by comparing our atomistic tight binding re- 
sults with room temperature photo-modulated re- 
flectance measurements. Our theoretical results are 
in good agreement with the experimental data and re- 
produce the crossover at 9% Bi composition, beyond 
which Aso exceeds the energy gap E g . Our theoreti- 
cal results explain that for the higher Bi composition 
samples with Bi fractions of 8.5% and 10.4%, the pres- 
ence of disorder, including a large number of Bi pair 
and cluster states close to the host matrix light-hole 
band edge, signifcantly increases the mixing of the 
LH band edge with Bi-related resonant states, and 
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therefore makes it difficult to identify a well defined 
LH band edge energy at higher compositions in the 
strained alloy. This is consistent with the poor fit ob- 
tained to the PR spectra above the band edge at these 
compositions. We have also analysed the impact of su- 
percell size and of statistical variations in the Bi distri- 
bution on the calculated electronic structure, showing 
that large (> 4000 atom) supercells are required to 
describe the alloy band structure, and to understand 
the experimentally observed inhomogeneous broad- 
ening of the transition energies. By comparing the 
band edge shifts in compressively strained supercells 
with those calculated for free-standing (unstrained) 
supercells, we have derived composition dependent de- 
formation potentials for the GaBi^Asi-a alloy band 
edges. Finally, our calculations suggest a type I band 
alignment at the interface between GaBi^Asi-z and 
GaAs, in accordance with recent experimental anal- 
ysis. Overall the calculations are in good agreement 
with the PR measurements and provide useful new in- 
sight for the development of this interesting new class 



of semiconductor materials. 
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